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Abstract 

The nonlinear dynamics of thermal and electromagnetic perturbations in the vortex 
state of type II superconductors is analyzed with account of dissipation and disper- 
sion effects. A theoretical analysis shows that nonlinear thermal and electromagnetic 
dissipative waves (structures) may be formed under certain conditions on the sam- 
ple surface. The structures possess a finite-amplitude and propagate at a constant 
velocity. The appearance of these structures is qualitatively described and the wave 
propagation velocity is estimated. The possibility of experimental observation of 
thermomagnetic shock waves is briefly discussed. 



Introduction 



Recent years a greater number of comparatively examples of spontaneous origin of spa- 
tial and temporary dissipative waves (structures) in unordered systems - self-organization 
processes - become well known. By the terminology of Prigozhin [1] dissipative structures 
are organized condition in the space and time, which can move over to the thermodynamic 
balance condition of only by jumping (as a result of kinetic phase transition). If the de- 
flection of a nonlinear system from the balance exceeds a critical value, these conditions 
can become unstable. In this case the system moves over to the new mode and becomes 
a dissipative structure, which appears and exists due to the dissipative processes. There 
are different types of dissipative structures [2, 3]. 

At present, superconducting systems with high critical field strengths and current 
densities are widely implemented in various advanced technologies. However, successful 
operation of the superconducting materials is only possible provided that special measures 
are taken to prevent a system from the thermal or magnetic breakage of superconductivity 
and the transition to a normal state. For this reason, one of the main problems in the 
investigation of properties of superconductors is predicting of the superconducting state 
breakage caused by dissipative and nonlinear effects related to viscous motions of the 
magnetic flux. This explains the considerable interest in the study of dissipative and 
nonlinear effects in superconductors that has arisen in recent years. 

The existence of essential nonlinear and dissipative effects, connected with Joule heat- 
ing at the moving of magnetic flow, creates different regimes of automodel processes, 
describing the evolution of thermal and electromagnetic perturbations in the vortex (re- 
sistive) state of type-II superconductor. One of measured regimes is the propagation of 
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thermomagnetic waves (normal region - the region heated with the temperature, which 
is higher than the critical temperature Tc) [3]. In early investigations [3, 4] it was shown 
that the nonlinear stage of evolution of thermal and electromagnetic perturbations in 
superconductors is determined by the stationary profile in the form of a running wave. 

In the present work the nonlinear dynamics of thermal and electromagnetic pertur- 
bations evolution in the vortex state of type II superconductors is investigated with the 
account of dissipation and dispersion effects. It was found the existence of nonlinear ther- 
momagnetic waves, describing the final stage of evolution of thermal and electromagnetic 
perturbations in the vortex state of superconductors. An estimate for the width and the 
velocity of a propagating wave is obtained. 

§1. Running wave in the vortex state of superconductors 



The evolution of the thermal (T) and electromagnetic (E, H) perturbations in the 
vortex state of superconductors is described by the nonlinear equation of the thermal 
conductivity [4] 

i.— ^W[KWT]+jE, (1) 

by the Maxwell equations 

r-o(/f=|j (3) 

and by the equation of the vortex state 

j^l{T,H)+l{E), (4) 

here u = i^iT) and k = k{T) are the heat capacity and the thermal conductivity respec- 
tively; jc = jc{T, H) is the critical current density and > = jr{E) is the resistive current 
density. 

The above system is essentially nonlinear because the right-hand part of Eq.(l) con- 
tains a term describing the Joule heating in the region of the resistive phase. Such a 
set (l)-(4) of nonlinear parabolic differential equations in partial derivatives has no exact 
analytical solution. 

Let us consider a planar semi-infinite sample {x > 0) placed in external magnetic field 
-> dH 

H — (0, 0, He) growing at a constant rate ^— =const. According to the Maxwell equation 

at 

(2), there is a vortex electric field E — {0,Ee,0) in the sample, directed parallel to the 
current density j: E \\ j; here He is the amplitude of the external magnetic field and Ee 
is the amplitude of the external electric field. 

In this section we will discuss the question on possibility of existence of automodel 
solution to the initial system of equations (1) - (4) in private derivatives with independent 
variables x and t in the form of running impulse, describing the nonlinear evolution of 
thermal and electromagnetic perturbations in superconductors. 
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The system of equations, describing the dynamics of evolution of thermomagnetic 
perturbations in the vortex state of superconductors is essentially nonlinear and it is 
impossible to integrate them in this form. As well known that, it is difficult to solve the 
problem with initial and boundary conditions for nonlinear differential equations in private 
derivatives by common methods, so that we can not use the superposition principle. 
Therefore it is reasonable to search other ways of the solution. There is such method, 
which allows the transformation of resemblance, due to this the system of equations 
with private derivatives can be reduced to the system of common differential equations. 
Recent years more attention is given to investigation of such automodel solution [5] . The 
obtained by the method differential equations with common derivatives can be integrated 
in the close form or reduced to investigation of one equation, which can be solved either 
analytically or numerically. 

In many problems of mathematical physics the automodel solution of running wave 
type is of great interest [5]. For important class of stationary running waves the distri- 
bution of characteristics of waves is still unchangeable on time. In accordance with this 
determination the solution in the form of running wave can be presented in the form 



where u - the characteristics of the event under consideration, x - coordinate, t - time, v 
- the velocity of the running wave. The solution (5) presents a wave, running either in 
negative or positive direction of x with constant velocity v. 

The possibility of reducing of the system of differential equations in private derivatives 
to the system of common differential equations in automodel derivatives greatly simplifies 
the problem in the mathematical point of view and in many cases allows to find exact 
analytical solution. Let us consider the problem about in what conditions the system of 
initial differential equations in private derivatives with independent variables x and t can 
be reduced to the system of differential equations in common derivatives, depending on 
only one variable ^(x, t). The solution of the system of equations (l)-(4) may be presented 
as a function of new automodel variable ^(x, t): 



Substituting (6) into the system (l)-(4) gives, as a result of simple differentiation, the 
following system 



u = u{x ± vt), 



(5) 



T = T[i{x,t)i 



E^E[i{x,t)i 



(6) 




dx^ d^ ^ \dx J de c2 [dT di ^ dE d^ dt' 



(7) 



(8) 
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In order the system (7), (8) was only function from ^ at the substation (6) it is required 
carrying out the following conditions: 

where A,C,G are functions to be determined below. Solving first two system of equations 
(9)-(ll) we have a relation 

G«)f^^(0f. (12) 
Hence just follows relationship between G and A 

GiO = ^^(0 , (13) 

here -u is a free constant of integrating of the equation (12). From (9) and (10) follows 
that i{x,t) must satisfy single-line equation in private derivation 

di .... 
— — u— (14) 
dt dx ^ ^ 

the only solution of which is the function 

^{x,t) = F{x-ut). (15) 
Using (15) we can immediately obtain 

G(0 = 1, A{0^-Fu, C(0 = 0. (16) 

It is possible to ensure F = 1 by the transformation of coordinates and times. Thereby, 
we find final automodel substitution 

C^x-ut, (17) 
corresponding to solution of a running type wave [5] . 

§2. Differential equations for nonlinear E and H waves. 



For the automodeling solution of the form (17), describing a running wave moving at 
a constant velocity v along the x axis, the system of equations (l)-(4) takes the following 
form 



dT d 



+ be(T)+J,(E)]E, 



(18) 
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(PE 

de ^ 

dE_ 



4:71V dj 
vdH 



(19) 
(20) 



The thermal and electrodynamics boundary conditions for equations (18)- (20) are as 
follows: 

T{i ^ +oo) = To, — (e ^ -oo) = 0, 



(21) 



-oo 



£;(e^+oo) = o, E{^- 

here Tq is the temperature of the cooling medium. 

In this section in view of considerable analytical difficulties involved in solving the 
exact problem, we will restrict our consideration to the so-called Bean-London critical 
state model [6] by assuming that j — jc{T, He) — jc{T) and 



jciT) = Jo- a{T -To), 



(22) 



where jo is the equilibrium current density, a is the thermal heat softening coefficient of 
the magnetic flux pinning force [7]. 

The characteristic field dependence of jr{E) in the region of sufficiently strong electric 
field {E > Ef) can be approximated by a piece-wise linear function jj. pa afE, where 



H. 



C2 



H 



is the effective conductivity in the flux-flow regime; 77 is the vis- 
7r/ic 

cous coefficient, $0 = is the magnetic flux quantum, (T„ is the conductivity in the 

2e 

normal state; is the upper critical magnetic field, Ef is the boundary of the linear 
area in the voltage-current characteristics of the sample [8] . In the region of small field 

E 

strengths {E < Ef), we assume a relationship jr{E) fa ji In — to be vahd with ji being a 

Eq 

characteristic local current density scatter (related to a pinning force inhomogeneity [9] ) 
on the order of ji ~ 0, Oljc; -E'o=const. This relationship between j and E is due to a 

thermoactivated fiux-creep. The thermoactivatcd flux-creep is a principal mechanism of 
the energy dissipation during the magnetic flux penetrates deeply into the sample. 

Exchiding variables T(^) and H{C,) from Eqs. (18) and (20) , and taking into account 
the boundary conditions (21), we obtain an equation describing the electric field E{^) 
distribution (£^-wave): 



d'^E 



de 



+ 



Attv djr dE ATTv^a N{T) - N{To) 



dE d^ 



k{T) 



aE^ 
2k{T) 



= 0, 



(23) 



dE^ 



where the dependency T = T ^E, j is defined by expression (18), (19) and have the 
form 



jo+jr{E) + 



Anv dE 
~^~d^ 



(24) 
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Here N{T) = Ju{T)dT. The velocity of E-wave is determined by the Eq. (23) with 



account of the boundary conditions 



El 



Stt 



N 



ro + -[je(r)+j,(^)] 



iv(ro) 



(25) 



Using Eqs. (20) and (23) one can easily find the expression for distribution of the magnetic 
field H in the case of H-wme 



(PH 4nv 



djr 

dE 



dH N(T) - N(To] 
+ ca- 



av H"^ 



0. 



(26) 



The velocity vh of H-wave is connected with its amplitude Hg by the following expression 

(27) 



N{T)-N{To) = ^. 



The expression (27) shows the adiabatic nature of the H - wave propagation: the magnetic 
energy, accumulated during decay of the field, is spent for heating the superconductor in 
the region close to the wave front. Therefore, depending on the character of the external 
electrodynamics conditions on the surface of the sample there may exist two type (E and 
H) of nonlinear stationary waves in the superconductor. 



§3. Thermomagnetic shock waves 



The first integrals of equations (18)- (20) have the form 

c2 dE 

E=-H, (30) 

c 

(an integrating constants in Eqs. (28)-(30) are equal to zero on the strength of the 
boundary conditions at ^ = —oo). Using integrals (28), (29) and excluding variable T{^) 
we can get the following differential equation for the distribution of E{^): 

^,,a,.)f,,v.-|-^o. (31) 

Here the following dimensionless parameters are introduced: z — —, /3 = — ^, L — 

H D 

is the depth of the magnetic field penetration into the superconductor, r = — - is 

47rjo D. 



m 



6 



the parameter describing the ratio of the thermal = — to the magnetic 

diffusion coefficient, = is the time of thermal diffusion, and E^^ = —— is a constant 

parameter. 

One can see that the mathematical model described by Eq. (31) can be considered as 
a decreasing nonlinear oscillator at the presence of a friction force 

F.-P(^^r)§, (32) 

(where z is the analogue of time, E is the coordinate of a " material point" ) with a potential 
[5] 

7712 



The analysis of the phase plane of Eq. (31) shows that there are two equilibrium points 
in it: Eq is a stable node and Ei is a saddle node. These two equilibrium points (£"0, -E-i) 
are separated in the phase plane by the separatrix AB (Fig.l). At z ^ +00 the material 
point is in the point Ei, and at 2; — ^ — cxd it moves to the point Eq. The transition from 
one equilibrium state to another one occurs only monotonycally. This solution by joining 
these two equilibrium points may be presented in the form of the shock-wave-type with 
amplitude Eg running at a constant velocity ve (Fig. 2). 

Let us consider the solution of the Eq (31), corresponding to the limiting case r << 1 
(hard superconductors). It follows from the r << 1 condition that the thermomagnetic 
perturbation develops nearly adiabatically [8]. Therefore in this approximation the ther- 
mal conductivity is negligible and corresponding terms, containing dissipative effects can 
be neglected in Eqs. (28)-(30). Then for the linear region of the current- voltage charac- 
teristic the solution to Eq. (31) can be presented in the form 



1 , 
tanh— -Iz — zn 

2 1+rr ^ 



(34) 



Here Zq is a constant. The condition r << 1 means that the magnetic flux diffusion 

Di is considerable faster than that of the thermal diffusion D^. Then we can assume 

that the spatial scale of the variation of the magnetic field penetration is essentially 

(X^E 

larger than the corresponding thermal scale Lt, therefore the spatial derivatives — — 

iXZ 

contain the small parameter -— << 1. It is easy to check fairness of this approach by 

Le 

differentiating Eq. (31): 

(fE(dEy^ 

r 

Maximum inaccuracy of given approximations is of the order of -—. For example, 

(1 + tY 

for T = 1 it is 25 %, and in the limiting case r — > or r — > 00 is very small. 

The expression (34) describes the structure of the thermomagnetic shock wave, pene- 
trating into the superconductor (see. Fig. 2). Using the boundary conditions for E-w&ye 
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one can easily find the wave velocity 



L 

VE=- 



2tE^ 



1/2 



with amplitude E^ and with front width 



5z= 16 



1 + T 
rV2 



E,, 



1/2 



(36) 



(37) 



sm 



Numerical estimations give the value ve = 10^ -i- 10"^ • — and 5^ = 10 ^-^10 ^for 

sec 

T = 1. 

Thus, the nonlinear stage of evolution of thermal and electromagnetic perturbations 
in the vortex state of superconductors is determined by the stationary profile in the 
form of running wave. It seen that, the transition from superconducting state to normal 
can occur by the expanding of stationary thcrmomagnctic wave, which is caused by the 
balance between nonlinear and dissipative effects, caused by Joule heating at the viscous 
motion of Abrikosov vortices inside superconductor. 



§4. NonlineEir waves in the flux-creep regime of heird superconductors 



As mentioned in Sec. 2 in the viscous fiux-fiow regime the electric conductivity cx/ is 
independent of the electric field the profile of thermomagnetic waves is independent of 
the background electric field in the sample. In the region of small electric field strengths 
the current-voltage characteristic of the sample has a nonlinear part and it is determined 
by the following logarithmic dependence 



3r{E) 



■ 1 ^ 

£/0 



The current-voltage characteristics of hard superconductor is described by the last equa- 
tion over a wide interval of temperatures T and magnetic fields H. This factor may 
considerably affect the character of nonlinear thermomagnetic wave propagation in the 
sample. In this section we consider a problem of thermomagnetic wave profile, with ac- 
count the nonlinear current-voltage curve of superconductor. Taking into account the 
nonlinear dependence jr{E), we can represented Eq. (24) in the following form 



a 



, E 47rvdE 
Jo + Jiln— + ^-7- 
Eo c^ dz 



Substituting the last expression into Eq. (18), and taking into account the boundary 
conditions (21), we obtain an equation describing the distribution of - wave: 



d^E 



+ /9 



1 + 



<Ta{E)E 



dE 
dz 



+ (3't 



Jo 



ME) 



+ 



ad{E) Eq 



E^ 
2K 



Here = (Jii{E) is the differential conductivity in the flux-creep regime 
sponding equation of state is obtained using the relationship [10] 



flE^ = X{E) 



1 + — In — , 
Jo Eq 



(38) 
The corre- 

(39) 
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where fl — 2jocJ {Ee)EK(3'^T is a constant. As seen from the plots of X[E) versus E 
(Fig.3), there exists a single intersection point of the curves y — QE'^ and y — 1 + — In — , 

Jo Eq 

which corresponds to a single stable equihbrium state. The stability of this state is 

(PE 

determined by the sign of the derivative ——- in the vicinity of the equilibrium point. The 

wave velocity ve can be determined from Eq. (39) using for the boundary conditions and 
have the form of 



Ve 



f 



2r 



joEn 



-1/2 



1 + ^ In . (40) 

Now we can use Eq. (30) and readily derive an equation describing the field distribution 
for a nonlinear if-wave. The wave velocity Vh of the H-wave is given by the formula 



cE 

vh^ — exp 



ad{E) 
2tji 



E. 



Ctn 



Jo 



(41) 



It is seen that the wave velocity increases exponentially with the amplitude Hg. For a 
sufficiently small amplitude H^ < Ha, where 



Hn 



L 



'iTjoEn 

<Jd{E) 



1/2 



(42) 



the wave velocity is negligible small, which corresponds to the case of the thermoactivated 
flux creep. For Hg — Ha the wave propagates at a finite constant velocity vh- The 
expression 



E — Eq exp 



aa{E)^ (LH, 



-E, 



Jo 
ji 



(43) 



2Tji " V ct^ 

represents the intensity of the spontaneous electrical field caused by the thermally acti- 
vated magnetic fiux-creep inside the region warmed up during wave propagation. In this 
case the maximum heating of the superconductor in the region immediately close to the 
wave front is described by the relationship 



T-Tq 
To 



e 



H^. 



Stti/To 



Let us now estimate the maximum heating from the following typical values of the 
physical parameters of the sample; iJg ~ 10^ • Gs, v — 10^ • erg/(cm^ • K), Tq = 4.5 • K, 
we have © fa 0.2 ~ 0.5. At this value of the heating the wave velocity is of the order 
of I'ij ~ 10 ~ 100 ■ sm/sek. Finally it should be note that taking into account the 
nonlinear dependence of the current density j on the electric field strength E should not 
qualitatively change the main results, since the character of the equilibrium state on the 
phase plane remains essentially the same. 



§5. The structure of waves 



The nature of the structure of the nonlinear E (or H) - wave can be determined 
by investigating the asymptotical behaviour at z — > ±oo of solutions for E{z) in the 
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close vicinity of the equilibrium points E = E^. Linearizing the equation (38) near the 
equilibrium point (weak-nonlinear wave) it can be represented in the following form 



(fE ^dE 



dE^ 



<yd{Ee)E, 
'2jo{To,E,) 



(E E,) 
2E. 



0. 



(44) 



The analysis of the phase plane (E, ) of Eq.(44) shows that there is a single equilibrium 

dz 

point, which is a saddle under the condition 



Jo(To,^e) > lfd{E,)Ee. 

The spatial scale of the variation of the solutions of Eq.(44), obviously, determines the 
width of the wave front A^o^ 



5zq 



1 - 



<yd{Ee)E, 

2io(To,£;e) 



n -1 



(45) 



An asymptotical solution of the Eq. (44) can be found in the field of 2; — > -|-oo(£' — > 0) 
using the voltage-current characteristics 



J ^ E^'^- {n » 1), 
in the flux-creep regime E « Eq ol the following form 

E{z) = E^e-xjp{-nl3z). 



(46) 



(47) 



The total width of the wave front is determined obviously, by the value A^^o since in the 
region of 2; > (on the strength of the condition n >> 1) the amphtude of the wave 



sharply tends to zero (on the scale Az* 



nl3 



« A.zo). It is seen that the width of the 



E-wecve front is determined mainly by the largest of the parameter heat capacity v and 
thermal conductivity k coefficients (/3 ~ — ). 



§6. Waves with the finite amplitude 



It is notice able that the evolution of perturbations of temperature T{x, t) and fields 

E{x,t) and H{x,t) essentially is determined by the equation of resistive state (4) and 

physical parameters of the sample. It should be noted that inclusion of the temperature 

dependences of the parameters k and v substantially complicates analytical calculations 

of the wave propagation dynamics that is described by the system of Eqs. (28)- (30). 

In most cases, the changes in the local values of these parameters in the sample can be 

considered negligible comparatively to the characteristic scale of temperature variations. 

Hence, we can take these parameters to be constant. Indeed, the investigation revealed 

that the thermal conductivity almost does not affect the character of the stationary wave 

dT 

propagation. This stems from the fact that the thermal flux k— vanishes at stationary 
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points of the system a.t z ^ ±00. However, the temperature dependence of the heat 

■ ^ T 3 

capacity should be taken into account. Such a dependence is represented as ly ~ uq 



To 

over a wide range of temperatures [8]. Then Eq. (28) can be written in the following form 



-i/ov 



T-TqA^ dT -2 



4Tn 



dz Snv 



By eliminating variable T{z) from the last relationship and employing the boundary 
conditions (21), we obtain 



dz"^ 



27r 



En & 



(TfE 



dE' 



aTo inavTo dz 



1 



^ dE E^ 
+ f3T— = 



dz 2E, 



(48) 



According to the qualitative theory [11], the equihbrium states are found from the condi- 
tion 



1 + 



OfEV 
aTc) J 



- 1 



E' 



(49) 



An evident property of system (49) is the absence of closed curves that are fully composed 

of the phase trajectories in the phase plane ^E, -j— = P^. The proof of this statement 

can be based on the Bendixson criterion [10]. The number of stationary points (one or 
three) and their type are determined by the parameter 



W = 2TTiyoTo- 



(50) 



E2 correspond to the condition 
at W > Wk. The 



The three equilibrium points E = 0, E = Ei and E 
W < Wk = - (Fig. 4). There is only one singular point Eo 
parabola and the quadratic curve in Eq. (49) are tangent at = Wk, i.e., this condition 
corresponds to the coincidence Ei — E2 — E* ° 



yields the following waves: 

1) Ei,2 = E*[l + 2,2{Wk-Wy/^] 



7 af 



at 



The direct solution of Eq. (49) 



Wk-W 



« 1 



2) El = 87r 



(51) 



Eo 



at Wk»W. 



Analysis of the phase plane {E, P) shows that the points Eq = and E = E2 are stable 
nodes and that £^ = £"1 is a saddle. In addition to the separatrix EiEo, the set (49) has 
the separatrix E1E2 connecting the points Ei and E2E1. This means that two types of 
waves with amphtudes AE — Ei and AE — E2 — Ei can exist in the superconductor 
(Fig.5). Evidently, wave I has an amphtude of the order E^ ai W ^ W^, its velocity 
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is determined by equality (49) aX E = Ey. Equation (49) has two stationary points at 
W « Wk'- £^0 = is a stable node and Ei = 2(3^tE^ is a saddle. The separatrix that 
connects these two equihbrium states corresponds to a "difference" -type solution with 
amplitude E^., which is related to the wave velocity ve and the wave front width Az by 
the equations (36) and (37), respectively. Wave II has a small amplitude at 1^ ^ Wk 

-—^AA{Wk-Wy/^«l, (52) 

and its velocity is inversely proportional to the amplitude at W « Wk- Such an exotic 
dependence oi ve on AE — Eg most likely means that the waves of this type are unsta- 
ble. Finally it is notice able that observation of the second-type waves becomes possible 

in finite-sized samples with asymmetric boundary conditions. The above investigations 
prove the possibility of applying the results obtained to high-temperature superconduc- 
tors cooled to liquid-nitrogen temperatures (T = 77 K), providing that the values of the 
physical parameters of the sample are known. 

§7. Stability of nonlineEir shock waves 



In Section 3, it was demonstrated that, depending on the surface conditions, stationary 
nonlinear thermomagnetic shock waves of two types, E and H, may exist in the supercon- 
ductor sample. In this Section the stability of nonlinear shock waves with respect to small 
thermal and electromagnetic perturbations in a hard and composite superconductors is 
studied. It is shown that spatially bounded solutions may correspond only to the per- 
turbations decaying with time, which implies stability of the nonlinear thermomagnetic 
wave. 

In order to study the stability of a nonlinear wave with respect to small perturbations, 
it is convenient to write a solution to Eqs. (l)-(4) in the following form: 



T{z, t) = T{z) + ST{z, t) exp 
E{z, t) = E{z) + SE{z, t) exp 



At 



(53) 



Here T{z) and E{z) are the stationary solutions and 6T, 6E are small perturbations; A is a 
parameter to be determined. Prom solution (53), one can see that the characteristic time of 
thermal and electromagnetic small perturbations tj is of the order of t^/A. Substituting 

(ST SE \ 
— r^, — TT << 1 1 
T{x) E[x) J 

obtain a differential equations for the distributions 6T and 6E in the following form 



v^JT =^^+[i(z)+ a,E(z)] SE - aE(z)ST , 



L2 cH^ 



(54) 
(55) 
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Eliminating the variable 5T by using the relationship (55) and substituting into Eq.(54), 
we obtain a fourth-order differential equation with variable coefficients for the distribution 
of small electromagnetic perturbation 5E: 



~d?~^^~d^ 



E, 



K 



^ - + A[At - B{z)]SE = 0, (56) 



Here B{z) = — — r characterizes the spontaneous heating of superconductor caused by 

a small perturbations of 6T and 6E. The condition of existence of a non-trivial solutions 
of Eq. (56) combined with boundary conditions (21) allows to define the spectrum of 
eigenvalues of A and the stability of the thermomgnetic wave in the sample, accordingly. 
From the expression (56) just follows that if ReA < 0, then spatially bounded solutions 
of 6T and 6E at infinity corresponds only a damping perturbations on the time. This 
problem is complicated, and its analytical solution cannot be found in a closed form. Let 
us consider the problem of the thermomagnetic wave stability in the limiting case of the 
hard (r << 1) and the composite (r >> 1) superconductors, which is presented a greatest 
practical interest. 

§7a. Hard superconductors (r << 1) 



As well known (sec. e.g., [8]), the magnetic fiux variations in a hard superconductors 
occur at a much greater rate as compared to those of the heat transfer, so that r << 1. 
The adiabatic character of the perturbations development leads to the predominance of 
magnetic fiux diffusion over heat diffusion in the sample, Dt « Dm- In this 
seen from solution (53), the characteristic times of temperature and electromagnetic 
field tm perturbations have to satisfy the inequalities tj << (A >> 1) and tj » tm 



(At << 1). Assuming that r << 1 in the limit of A >> , which corresponds to a 

" fast" instability [4] , we obtain the equation for determining the eigenvalues of A 

^ + /5-^ + [BW-Arl5B = 0, (57) 

In hard superconductors with r << 1, the instability threshold depends on the ther- 
mal boundary conditions at the surface of the sample only slightly. Therefore, the 
electrodynamic boundary conditions at the boundaries of the current-carrying layer 
{z = +00, z = —00) can be neglected and one can keep only the thermal boundary 
conditions to Eq. (57). Using the substitution of variables of the type 

dE _ iI-tE,^ 



Eq. (57) can be rewritten in the following form [12] 



di2 + 



0. (59) 
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Here 



^(2;) = SEexp 



y ^ —{z- Zo), A = At. 
The Eq.(59) with the help substitution 

e 

■iP{x) = {1- x'^)2^{xy, l-x^2s, 

X — thy, 1 — = ch~^y, 
e = VA, P = 1, 
can be reduced to a standard hypergeometric equation [13] 

s{l-s)^-[e + l- 2s(e + 1)]^ - [e(e - 1) - p(p - 1)]$ = 0. (60) 

Here, the even and odd integrals can be presented as 

$1 =F(e-l,e + 2,e+l,s), (61) 

$2 = s-^F(-l,2,l-e,5), (62) 
where F is the hypergeometric function [13]. Returning to the variable y we have 

*i(y) = (chy)-F f e - 1, e + 2, e + 1, ^-^] (63) 



*2(y) = (thy - e) exp(ey)F |^-1, 2, 1 - e. 



thy^ 



(64) 



The function must be finite at the singular point s = 1. The e values for which $1 
is finite evidently correspond to a discrete spectrum: = 0, 1 or Aj = —t~^,0. The 
function $2 is finite only for e = 0. Any solution in the form of a running wave is 
characterized by translation symmetry, which implies that a "perturbed" stationary profile 
E{z) determined by Eq. (57) corresponds to the eigenvalue of the ground state Aq = 0. 

dE 

Differentiating Eq. (31) with respect to z, one may readily see that — ; — is an eigenfunction 



corresponding to Aq = 0. Indeed, the perturbation 5E 



dEp 

dz 



dz 

essentially represents a small 

dEn 



wave displacement. Thus, we may suggest that tha function — — exponentially tends to 

dz 

zero for z +00 and the corresponding eigenvalue is zero. Therefore, the problem cannot 
possess positive eigenvalues and i?eAj < 0. This result implies that the wave is stable with 
respect to relatively small thermal 5T and electromagnetic 5E fluctuations. The analysis 
of the second linearly-independent solution (64) leads to the same conclusion. 



§7b. Composite superconductors (r >> 1) 
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In composite superconductors, an efficient value of effective electric conductivity aj is 
much above than in hard superconductors. As a result any intensity of the thermal and the 
electromagnetic small perturbations is valid an inequality r >> 1. It may be assumed that 
induced the normal current afE caused by raising of temperature compensates a droping 
of the critical current jc{T) and obviously, impediments of an entry of a magnetic flow into 
sample. On the other hand in composites the thermomagnetic instability development is 
accompanied by "slow" perturbations with characteristic time of a growth << tj << tm 

(A << 1) and consequently, >> t— = -f {^'t >> 1) [§]■ In the approximation A << 

At a 

from Eq. (56) we obtain the following expressions for determining the eigenvalues A 

K 



(f6E 



+ 



/5- 



E{z) 



d6E 

dz 



A 1 + 



E{zy 



SE^O, 



(65) 



The last equation with the help substitution 



^{z)^SE{z)chj{z-Zo), 



can be reduced to a canonical form [12] 



[H - A]* = 0, 
U{y)^-l + -^ + A,thj; y=^(z 



(66) 



zo), Ai = Xp- 



(67) 



The stationary wave is stable if the eigenvalues of the operator H include the negative 
Ai < 0. As a boundary conditions for Eq. (66) serves the condition that the solutions 
^(y) are bounded at infinity. The odd term in Eq. (67) describes the effect of the thermal 
mode on the dynamics of electromagnetic perturbations in the superconductor. Using the 
standard procedure [12] we can readily find the exact solution of Eq. (66) in the form 



1 



1 



^{y) = (1 - thy)""" 2(1 + thy)""" 2 F ( a, /5, 7, 



(68) 



Here a constant parameters a, 7 are determined by the following relationships 

a ^ p + q + 3, 

I3 = p + q, 
7 = 2p+ 1. 



and 



VI + 3Ai 



2 ' " 2 

It can be seen that the spectrum of eigenfunctions in (68) is continuous [12]. Analysis 
of the asymptotic behaviour of solution (68) at y — > ±00 shows that SE{y) is bounded 



15 



only at Ai < 0. Indeed, one of the partial solutions of Eq. (66), remaining bounded at 
y ^ +00, has the form 

5E(y) ^expiupv) X F{p + q + 2,p + q-l,2p+l,ex.p{-2y)); 1 - 2p ^ Up. (69) 

It is obvious that ai y ^ +00 (exp(— 2|/) — > 0) to the asymptotic (68) corresponds to 
exponentially decaying at positive cUp and growing wave as exp^iujpy) at negative cUp. It 
may be assumed that the obtained solution exhibits regular asymptotic behaviour at 
y +CXO i.. the function 5E{y) is bounded only at those eigenvalues for which Ai < 0. 
Let us represent the second linear-independent solution of Eq. (66) as: 

F(p + g + 2,p + g - 1, 2p + 1, s) = ^-^^^^^it^^L^x 

[p- q-l){p- q + 2) 

xF(p + g + 2,p + g-l,2p + l, l-s)+ (70) 

+ r - -q-hp-q + 2,-2q + l,l-s). 

{p + q + 2){p + q + l) 

where is the Euler's function [13]. At y — > —00 the following asymptotic representation 
is valid 

5E{y) ^ exp(t^qy); 1 - 2q = cUq. (71) 

It can be seen that the space-limited thermomagnetic perturbations are damped both at 
(Vp > and at a;^ > i.e., Ai < 0. Let us consider a some specific cases which follows 
from the properties of the asymptotic solutions of Eq. (66). 

a) — < Ai < 0; then the parameters a;^ > and a;^ > can be chosen substantially. 

3 

According to the (53) small perturbations damped exponentially with the passage of a 
time (see, [4]) (Fig.6). 

b) — - > Ai > 0; in this specific case the wave amplitude at the sufficiently large positive 

values y is damped exponentially and at the sufficiently large negative values of y it has 

27r 

the oscillating-damping profile with the wavelength — (Fig. 7). 

Jm • q 

c) Ai < —1; in this case to the solution corresponds oscillating-damping wave at y — > ±00 

27r 2n 

with the wavelength Lg = and Lp = (Fig. 8). 

Jm • q Jm • p 

Thus it is shown that only damped perturbations correspond to space-limited solu- 
tions, which means that the nonlinear wave is stable. 

Conclusion 



In conclusion note that the existence of essential nonlinear and dissipative effects, 
connected with Joule heating at the moving of magnetic flow, creates different regimes 
of automodcl processes, describing the evolution of perturbations in the vortex state of 
type-II superconductor. One of measured regimes is the propagating of thermomagnetic 
waves (normal region - the region heated with the temperature, which is higher that the 
critical temperature Tc) . The nonlinear stage of evolution of thermal and electromagnetic 
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perturbations in superconductors is determined by the stationary profile in the form of 
a running wave. It is seen that the transition from the superconducting state to the 
normal state can occur by expansion of a stationary thermomagnetic wave, which is 
caused by the balance between nonlinear effects, caused by Joule heating and dissipative 
effects. The results obtained make it possible to describe the final stage of evolution 
of the thermomagnetic instability in the vortex state of typc-II superconductors. It is 
possible that two kinds of stationary thermomagnetic shock waves to exist, depending on 
the electrodynamic boundary conditions on the sample surface. 

We should emphasis that this result was obtained for an arbitrary temperature de- 
pendence of the thermophysical parameters v and k of superconducting material and for 
an arbitrary function j{H). Moreover, since the system of equations (l)-(4) is invariant 
with respect to an arbitrary translation, the wave propagation condition can be found for 
an arbitrary critical current density dependence on T and H. 

Let us discuss the validity of the present approach. The applicability of the macro- 
scopic description is ensured by the inequality » d, d being the average distance 
between the Abrikosov vortices. The boundary conditions used are valid to describe the 
situation for a sample with a finite thickness 2L provided A,^ « L. It is seen that the 
shock waves under consideration can be observed only for sufficiently thick samples. 

The result obtained enables us to propose a possible version of thermomagnetic insta- 
bility development in superconductors. The initial stage of flux jump is characterized by 
the exponential growth of both thermal and electromagnetic perturbations. The insta- 
bility threshold as well as the instability increment are defined on the basis of the linear 
theory [14]. At a further stage the thermomagnetic shock 'H-wave' with the amplitude 
He = Hj spreads inside the superconductor. Such a process results in a total penetration 
of the magnetic flux into the sample. The velocity vh to be found from equation (41), 
determines the time of the wave motion inside the sample of thickness 2L 

Vh 

However, the normal transition occurring within the warmed region complicates the con- 
dition of shock-wave propagation. A more detailed description of a such two-phase con- 
figuration will be presented in a further report. 

Finally, let us briefiy discuss the possibility of experimental observation of thermo- 
magnetic shock waves. The necessary condition to observe the phenomena predicted is to 
ensure sudden initial magnetic flux lattice motion within the superconductor. For exam- 
ple, a thermomagnetic 'E-wave' could be generated by initiating the magnetic flux motion 

in the presence of an external magnetic field varying with the rate He- Then, a 'back- 

L ■ 

ground' electric field of order ~ —He arises inside the sample and directly determines 

c 

the amphtude of an 'E-wave' shock wave. We would point out that thermomagnetic non- 
linear waves may also be initiated by abrupt heat input leading to a sudden temperature 

increase ST = T* — Tq at the sample surface (accompanied by a simultaneous magnetic 
flux slip initiation inside the superconductor). The electrodynamic boundary conditions 
at the sample surface is determine the type {E or H - wave) of the propagating wave. 
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FIGURE CAPTIONS 



Fig.l. 


The phase plane of Eq. (31) . 


Fig 


.2. 


The structure of shock wave. 


Fig 


.3. 


The phase plane of Eq. (39). 


Fig 


.4. 


The equilibrium points of Eq. (49). 


Fig 


.5. 


The phase plane of Eq. (49). 


Fig 


.6. 


Asymptotical behaviour of waves at — < Ai < 0. 

3 


Fig 


.7. 


The wave amplitude at the sufficiently large positive values 



y is damped exponentially and at the sufficiently large negative 
values of y it has the oscillating- damping profile at — - > Ai > 0. 
Fig. 8. Asymptotical behaviour of waves at Ai < —1. 
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